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Abstract 

The effective actions of a scalar and massless spin-half field are de- 
termined as functions of the deformation of a symmetrically squashed 
three-sphere. The extreme oblate case is particularly examined as per- 
tinant to a high temperature statistical mechanical interpretation that 
may be relevant for the holographic principle. Interpreting the squash- 
ing parameter as a temperature, we find that the effective 'free ener- 
gies' on the three-sphere are mixtures of thermal two- sphere scalar s 
and spinors which, in the case of the spinor on the three-sphere, have 
the 'wrong' thermal periodicities. However the free energies do have 
the same leading high temperature forms as the standard free energies 
on the two-sphere. The next few terms in the high-temperature series 
are also evaluated and briefly compared with the Taub-Bolt-AdS bulk 
result. 
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1. Introduction 

The holographic 'principle' says, in its barest form, that the information con- 
tained in the interior of a space-time domain is encoded in a field theory residing on 
the boundary. In accordance with this, the gravitational entropies of the Taub-Nut- 
AdS and Taub-Bolt-AdS space-times have been computed by Hawking, Hunter and 
Page [1] and by Chamblin et al [2] with the aim of comparing with boundary (con- 
formal) field theories. The boundaries are symmetrically squashed three-spheres, 
possibly with identifications. 

The object of this paper is to detail the calculation of the effective action, i.e. 
the functional determinants, of scalar and spinor fields as functions of the squashing 
on the three-sphere. The effective action will be associated with a free energy and 
thence with an entropy. 

The fields will be free so this is only a prelude to a more realistic investigation 
and the present results may well have only a passing relevance to the holographic 
principle. Such a stopgap calculation has been suggested by Hawking et al in [1] . 
Nevertheless, the evaluation of the determinants is of interest in itself. Our results 
here will be confined to the determinants. Elaboration of the statistical mechanical 
interpretation is left for a later paper although some preliminary comparison with 
the bulk Taub-Bolt-AdS result is made. 

2. The basic situation 

The squashed 3-sphere appears as the spatial section of the frozen mixmaster 
universe and quantum field theory on this space-time was discussed by Hu, Fulling 
and Parker [3], Hu [4], Shen, Hu and O'Connor [5] and Critchley and Dowker [6]. It 
has been discussed in a Kaluza-Klein setting by Okada [7] and by Shen and Sobczyk 
[8]. 

A discussion of the vacuum energy of a massless spin-half field on the squashed 
3-sphere has been given in [9] and for this reason we will reexamine this case before 
looking at the scalar field. 

In reference [9], in order that rules of standard angular momentum theory 
should apply unmodified it was necessary to choose the radius of the unsquashed 
sphere, S 3 , to be a = 2. and the (standard) metric in Euler angles is 

ds 2 = (de 2 + sin 2 6d(j) 2 ) + l 2 (dt(j + cos6d(j)) 2 (1) 
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showing the (symmetrically) squashed S 3 as a twisted product, S 2 x S 1 . The circle 
has radius 2Z 3 , and hence a circumference of 4irl 3 . If, illustratively, this periodicity 
is translated into a temperature, we find (3 = 1/T = 47r/ 3 . Another way of saying 
this is to note that the range of ip is — > 4ir and then to interpret l 3 ifj as a Euclidean 
time. We will be particularly interested in the extreme oblate case, l 3 — > 0, when the 
metric reduces to that on the unit two-sphere. The relation between the (^-functions 
will constitute a type of twisted Selberg-Chowla formula. 

In the notation of Hawking, Hunter and Page [1], /| = E. It is possible to iden- 
tify points on the '^-circle', and this is the more interesting situation. Nevertheless 
it will not be considered here. The work [2] also treats only this simplest case. In 
the notation of [2], l\ — An 2 /I 2 . We denote the squashed sphere by S 3 . 

3. Spinor (— functions on S 3 . 

We ignore all restrictions on the boundary spinor theory arising from its em- 
bedding in some bulk theory. For odd-dimensional spaces the relevant operator is 
thus the Pauli one which, on the squashed three-sphere, is n = — za*Vj in terms of 
(covariant) Pauli matrices and the spinor covariant derivatives, which we will not 
exhibit here. The dimension of spinor space on S 3 is 2, which is the same as that 
on S 2 . 

The eigenvalues of n are determined to be (Hitchin [10], Gibbons [11], and [9]) 

"± = (2/s)- 1 Qll ± {n 2 + 4(/ 3 2 - l)q(n - q)) ^ (2) 

where the integers n and q emerge from the angular momentum quantum num- 
bers labelling the unperturbed states (before the necessary secular diagonalisation) . 
Thus n = 2L + 1 and q = n/2 — M with L an orbital label and M the projection of 
the total angular momentum. 

If Z 3 < 4, u + is positive and u_ negative, and for ease we will assume that this 
is so. The corresponding traced ^-functions are constructed separately, 

00 n 
n=l q=0 + 

and 

00 n—l 
n=2 q=l V -' 



2 



which exhibit the quantum number range restrictions with n the remaining degen- 
eracy of the modes. We note that, as a function of l 3 , nothing peculiar happens to 
the eigenvalues as Z 3 passes through 1 and also that there is a square root branch 
point at /| = 0. 

The (^-function for the squared operator IT 2 is 

Tt 3 ( 3 (s)=( + (2s)+(_(2s). 

The awkwardness of the eigenvalues, (2), restricted the analysis in [9] to a power 
series expansion in the squashing parameter l 3 . For the time being we shall continue 
with this expansion which, in any case, is adequate for the high temperature limit. 
The computation then reduces to that of the function, /(s), defined, for Res > 3/2, 
by 

oo n — 1 

/(^)=EE / 2 + fi*> where 5 2 = / 3 2 - 1. (3) 

In terms of f(s), Tr 3 ( 3 (s) reads, 
Tr 3 C 3 (s) =2(2/ 3 ) s (c(2s-l,/^ (4) 

This expansion is valid, numerically for l 3 < y/2. One can rejig the series to 
allow for larger values, but we will not bother. With this approach, it is not possible 
to discuss the l 3 — > oo limit. 

A preliminary aim is to calculate the (Euclidean) effective action, which we 
define to be Tr 3 C 3 (0)/2 = W sp . 

Noting that there is no conformal anomaly, Tr 3 £ 3 (0) = 0, the complete formal 
series is 



W sp = 2C'(-1, ll) - 2/ 2 C'(0, ll) + /'(0) + 4(2P + R) + -llf(2) + Oil?) (5) 
where, since f(s) has a pole at s = 1, we have to write 

f(s) = -^- + R + 0(s-l). 
s — 1 

It is expected that the interesting behaviour as Z 3 becomes small is contained 
in /'(0) and we must determine the analytic continuation of this quantity. Initially 
we are looking for terms which diverge as l 3 — > 0. 
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In [9] we employed the standard Plana summation formula and, for convenience, 
will proceed in the same way although there are other approaches. 

Extending again the q sum to to n in (3), an application of the Plana sum- 
mation formula yields, for deformations in the prolate direction, [9], 



f(s) = -C R (2s - 1) + ( R (2s - 2) / ^ ^ 

(6) 



r at r y> n 

' 7 exp(27rt) - 1 I ^ (n 2 + 45 2 (t 2 - itn)) 1 



(t - -f ) 



which is not a complete continuation as the n summation has yet to be dealt with. 
The pole at s = 1 is apparent giving P = — 1/2. 

For deformations in the oblate direction (0 < / 3 < 1), which we need for the 
high temperature series, singularities of the summand encroach into the relevant 
part of the complex q plane, i.e. < Keq < n. Concentrating now on this oblate 
case, 'extra' branch points occur at 

n , .n L n ,„ , . ^ 

q = — ±i ; = -(1 ±itan0) 

2 2 v^Tf 2^ ' 

and we give the details of the resulting contributions. Following the usual procedure 
(as in Lindelof [12]) we find (setting b 2 = — S 2 and omitting the n-summation for 
the moment) the additional pieces, 

dz n f dz n 



f dz n f 

Ju exp(-2niz) - 1 ( n 2 - Ab 2 z(n - z)Y J L 



u exp(-2niz) - 1 (n 2 - 4b 2 z(n - z)) s J L exp(2niz) - 1 (n 2 - 4b 2 z(n - z)) s 

where the infinite U contour runs anticlockwise around the cut from n tan^)/2 
to n/2 + zoo, while L runs clockwise around the corresponding cut in the lower half 
plane. Symmetry means that the two contributions are equal and we get, after a 
change of variable, z = (1 + i(ta,n9) n/2, the extra contribution to f(s), 



ita,n9 1 f d( 

If t[n 2s - 2 J c ((-l)"exp(n7rCtan^)-l)(l-C 2 ) s U 

where C is an anticlockwise contour running around the real £-axis, cut from 1 to 
oo. For s a non-negative integer, this expression vanishes, while, if s is a negative 
integer, it can be evaluated using residues, providing a useful numerical check. 
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Because of the exponential factor in the denominator, expression (7) converges 
for all s. It can be taken as the analytic continuation and should be added to (6), 
for deformations in the oblate direction, to give 



f( s ) =-C Ms- 1) +C??(2s-2) f - — 



2i 



DC 



dt f \ -v n 



-l\J2 ( n 2_ A h 2 (f 2 _ j±n) y (* "> £ ) f + (8) 



o exp(27rt) - 1 I ^ (n 2 - 46 2 (t 2 - itn))' 
itan6> 1 /" <iC 1 



*§' ^!^ 2s " 2 ic ((-l)"exp(mrCtan0)-l) (l-C 2 )" 

We shall denote the last two terms on the right-hand side of (8) by / 2 ( s ) an d 
/ 3 (s) respectively and consider f 3 first. 

If s is such that the integral converges at £ = 1, (7) can be converted to a real 
integral in the usual way by combining the upper and lower pieces of C. One gets, 
choosing appropriate phases, 

, / N n • / N tan6>>^ 1 f°° d( 1 . . 

fo(s) = 2 sin(7rs I— 75 — > — = / : : m (") 

If Ji ((-l) n exp(n 7 rCtan6') - 1) (C 2 - 1) 

In particular we can differentiate with respect to s to get the contribution to 
f(0), 



r°° dt 

^ Ji ((-l) n exp(n7rCtan^) - 1) 



(10) 



This is a simple example of a Lerch function. Expansion of the integrand allows 
the integration to be done to yield, from (5), the contribution to ^Tr 3 (' 3 (0), 



00 1 00 

W sp =2 J2 ^(-l) mn ne-™ ntane . (11) 

m=l n=l 

The sum over n produces 

w£Hp') = /a(o) = y — i y i r, (12) 

bp V ' 3K 1 £-{4ksmh 2 (k[3'/4) f^2(2A; + l)cosh 2 ((2A;+l)/3V8) V ' 

where, for convenience, we have denoted 47rtan6> by f3' . This expression is useful 
numerically. 
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For small Z 3 , the exponent in (11) can be written as — mn/3/4 with j3 defined 
by (3 = 47r/ 3 and so, dividing by (3 for normalisation purposes, we define the corre- 
sponding spinor 'free energy', 



^ 3 p ) (/3) = -^ip 3) (/3 / ) (13) 

We will postpone discussion of formula (13) until after the corresponding scalar 
case has been considered. 

However before proceeding to this, it is necessary to consider the other terms in 
(6) and (8) which are needed for the complete determination of the effective action 
as a function of the squashing. 

In [9], a further application of the Plana formula to the n-sum in (6) revealed 
a series of poles in f(s) at s = 3/2 — m, m = 1, 2, . . ., with residues 

r m = (- 1 ) m+l 22m "J r ( ( 7/2) 1/2) ( ^ " 1 ^ l " m ~ 2B ^ (14) 

where B 2m is a Bernoulli number (using the definition in Bateman [13] e.g. In [9] 
we used Lindelof's signs [12]). 

Here we will employ the Watson-Sommerfeld technique which is essentially 
equivalent to the Plana one. It has been used by Shen and Sobczyk [8] in the 
present context. 

Completing the square, we have, 

n 2 + 45 2 (t 2 -itn) = (n + iB) 2 + A 2 (15) 

where 

A 2 = ^(/g - l)t 2 = ~A 2 t 2 , B = 2(1 - l\)t = Bt with B 2 - A 2 = 2B. 



The signs etc. are appropriate for the prolate case. For oblate Z 3 < 1, we set 
A 2 = -C 2 . 



We leave the t-integration until last and just consider, for the oblate case first, 

En I I z 
r-Tf = — / dz- -irCOtllZ 
n=1 {{n + iBy-C 2 ) s 2zJ L (( z + iB) 2 -C 2 ) S (16) 

= - / <K , C ~ lB xs cot (tt(C - iB)) 
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where L is the anticlockwise contour surrounding the poles of the cotangent at 
C = n + iB, n = 1,2, — The £-plane has branch points at — C and +C. For 
symmetry's sake the associated cuts are arranged to run down the imaginary axis 
to (almost) the origin, and thence along the real axis to ±C. To make a choice, the 
+C-cut comes in from +ioo and the — C one from —zoo. The contour L is deformed 
to run down, just to the right of the imaginary axis, then around the right-hand 
side of the C-cut and finally to skirt the imaginary axis down to —zoo. We assume 
that s is such as to ensure convergence. 

Taking the phases into account, the integral along the imaginary £-axis is 



"2/ 



dy J+ciy e ~™ s cothn{y ~ B) ~IJ dy (yKc*y ^ a cothn{y + B) 



and remembering that, according to (8), we need 4 times the imaginary part of this 
expression we find, 



2 sin7rs 



/■oo rco y -\- B 

J dy (y2 + czy cothn{y ~ B) ~J dy ( y 2 + C 2y cothn{y + B > 



Now we make the split coth7rx = 1 + 2/(exp27rx — 1). The '1' part of this 
allows the integral to be done, 



2 sin7rs 



BC ^ns ~ 1/2) _ 2 r dy y + B 



F(s) J *(z, 2 + C 2 ) s (exp(27r(z/ + J B))-l) 

y-B 



+ 2/ dy- 



( y 2 + C ^y(exp(2n(y-B))-l), 

(17) 

and it is straightforward to check that the first term in (17) reproduces the poles 
in f(s) at s = 3/2 — to, with residues (14), after putting back the t-integration. 

The integral over the right-hand part of the C-cut can likewise be reduced to 
the form, if it converges, 

f° , 1 x sinh(27rS) - B sin(27rx) 
4sin7rs / dx— T —— — —-. (18) 



(C 2 - x 2 ) s cosh(27rS) - cos(2ttx) 
From these expressions we can derive /^(O), 

Ant 2 dt f f B . ... f° , xsmh(2nBt) - Bsm(2nxt) 



m = 



o exp(27rt) — 1 



' dyy coth(nyt) — / dx- 
Jo 



cosh.(2nBt) — cos(27rxt) 

(19) 
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We note that, for s a positive integer, the integral in (18) does not converge so that 
this method is not convenient for these values. 

The prolate case can be treated in the same way and gives 

,,, , AiT, , , f°° t 2 dt f^-z ydy . , 

In the near-round limit, l 3 ~ 1, /^(O) goes like (1 — /|)/3 in both cases. 

This approximation is most easily developed from the perturbation expansion 
of the original sum in (6). Since this provides a useful check, we give the first three 
terms obtained in this way. Defining l 3 = cosh 7 in the prolate and l 3 = coscj) in 
the oblate case, we find 

nf/ , 1 o /I 2tt 2 \ 4 / 2 4tt 2 16tt 4 \ fi . , 

and the corresponding oblate form obtained by setting 7 — > z0. Numerically, for 
7 = 0.5, (21) gives —0.123, in nice agreement with numerical integration applied to 
(20). 

To order l 3 , the complete effective action from (5) is 

W sp = - 2^C*(3) + /2(D) + &(0) + Oilt) (22) 

where /^(O) is given by (20) in the prolate and by (19) in the oblate case. In the 
oblate case, f 3 (0) is given by (11), otherwise it is zero. 

W sp , as given by (22), is the quantity we concentrate on. While the expres- 
sions in this section are numerically adequate, more suitable forms are developed in 
section 8. 

For future reference, we note that when A and C are zero, but B ^ 0, both 
(19) and (20) give the same result, namely just the first part of (19). One must, 
however, be cautious with the C — > limit in general, since the integrals in (17) 
may not converge. If one did wish to pursue this further, then a different contour 
deformation would be called for. Furthermore, it is possible to make a perturbation 
expansion in C, or in A, of the original summation in (6) or (8) the coefficients of 
which depend on the (imaginary part of) the ^function 

00 

En 
n (n + iB) s 

n=0 v ' 

similar to the Hurwitz (^-function with imaginary parameter. The general nature 
of this expansion does not distinguish between prolate and oblate deformations. 
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4. The scalar field 



We now turn to the massless scalar field. In the context of the squashed three- 
sphere, this has earlier been considered by Critchley and Dowker [6], Shen and 
Sobczyk [8] , and more recently by Shtykov and Vassilevich [14] who were concerned 
with the heat-kernel expansion for the Laplacian. We will use the scalar operator, 
—A + 1/4 for which the eigenvalues are 

X = ^q( n2 + 462{q + ^ {n ~ q ~ 1 2 ) ) (23) 

with degeneracy, n. The n label runs from 1 to oo, and q from to n — 1. 
We recall that the scalar curvature on S 3 is 

/ 2 

* = 2-| 

which correctly gives 3/2 when Z 3 = 1 for the round three-sphere (of radius 2) and 
2 for the two-sphere at Z 3 = 0. Our choice of operator restricts to —A + R/6 in 
the round case, and corresponds to a conformally invariant wave equation in four 
dimensions. It also reduces to — A + R/8 on the unit two-sphere when / 3 = 0. 

Other possible operators would be the minimal one, —A, and the one always 
conformally invariant in three dimensions —A + R/8. It is not difficult to accom- 
modate these by an expansion as in the spinor case, but since we are interested, for 
the moment, mostly in the leading high-temperature terms, the choice was dictated 
by the simplicity of the eigenvalues. 

The (^-function, J2^~ s = ( sc ( s )i * s 

oo n— 1 

CM = (2/ 3 ) 2s £ £ ^ (24) 

n=i g =o l n 2 + 46 2 (q+ \){n-q- \) J 

and the effective action is W sc = — Cg C (0)/2. 

The Plana summation formula is now applied, in the manner of [9] and dis- 
cussed also in [8] . We will consider the oblate case. Again extra singularities appear 
in the relevant band of the complex (/-plane. Replacing q by z and making the trans- 
formation of variable 

n — 1 .n 

z = h i — Ctant/ 

2 2 s 

the denominator becomes 

(iin 2 (i-c 2 )Y 
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showing the branch points at £ = ±1, as before. 

In this scalar case we use the alternative choice of summation contour, the 
band in the z plane being defined by —1/2 < Re z < (n — 1/2) (rather than to 
to — I). Then the full Plana formula, applied just to the q summation, is 



/cot 8 r 

2l{2h)2S So exp(27rt) + l { E ( n2 + 4^2 _ itn) y "ft-"*)}- 

z2 2s tan0V -J-^ / — 7 - 1 , s , 

^n 23 - 2 ^ ((-l)«exp(n7rCtan^) + l) (l - C 2 ) 

(25) 

Look at the last term first and again combine the upper and lower parts of the 
clockwise loop C to give 

-2 2s+1 sin(7rs)tan^y / y-— 7 - (26) 

V ; ^ x n 2s - 2 J x ((-l)-exp(n7rCtan^) + l) (C 2 -l) 

when the integration converges at Q = 1. 

The contribution to — (^ c (0)/2 from this term is then 



W s ^ = n tan 9 V n 2 / — - 

£i A (- 1 ) 



(27) 



L) n exp(mrCtan6>) + 1 

which can be written as 

Wi c 3 )(/3 / ) = jr ^ ^(_i)-(-+D ne - nm/3 '/ 4 , (28) 

and sums to 



m=l n=l 



^(3) / OA = * , i / 29 N 

S fct 8A: sinh 2 (fc/3V4) ^ 4(2fc + 1) cosh 2 ((2k + 1)0'/ 8) ' 

The corresponding contribution to the free energy <&H\[3) = ^Csc(O) is 

^c ) (/3) = -^s ( c 3) (/3 / ) (30) 

When l 3 is zero, the metric (1) reduces to that of the unit two-sphere. The 
scalar eigenvalues of —A + 1/4 on the unit two-sphere are n 2 /4 with odd n. The 
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degeneracy is n and so the standard finite temperature correction to the bosonic 
free energy reads 

m=l n=l,3,... 

The eigenvalues of the massless Dirac operator the unit S 2 are ±n, n = 1, . . . , oo 
with degeneracy 2n. So, (including positive and negative modes) 

^ m=l n=2,4,... 

The effective action densities on the squashed three-sphere, in the highly oblate 
direction seem to possess both boson and fermion characteristics. 
Split the sum in (28) into even and odd n, 



00 



1 00 00 ( -\\m 00 

Y- Y ne-^+yi^- Y ne"™^ 4 (33) 

m=l n=l,3,... m=l n=2,4,... 

so we can write, at least formally, 

= j (l^cWm - \K P wm)- (34) 

Hence the small (3 behaviour can be determined from known results on the two- 
sphere. 

The standard expressions for the Weyl terms, for thermal behaviour on a two- 
manifold, here the two-sphere, are 



TP' (R\ ~ Cb( 3 ) ^sc _ 2 Cr(3) 



and 

KM ~ -^P^ = -^F (36) 

where Cq c = \M\ and Cq P = |.M| Tr 1. The coefficients of the C 's and indeed all 
formal expressions, satisfy the functional relation 

Q' ap W) = Q' S cW)-2Q' sc m- (37) 
From (34) the leading terms are 

•2>Cfl~-^ + ^^ + ^ + ±ft(-i>- (38) 
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This result can be obtained more directly from the summed form (29), which 
differs from the spinor form, (12), just in the sign of the second term, and a factor of 
two. Moreover, this term has a logarithmic asymptotic dependence on f3 as f3 — > 0, 
which is expected since the C x heat-kernel coefficient on the two-sphere is non-zero. 
The corresponding sub- leading terms can be seen in the last two terms in (38) . The 
scaling of the In/? is provided by the size of the two-sphere (for which there is a 
conformal anomaly). The second term on the right-hand side of (38) is a simple 
consequence of the difference between 'effective' inverse temperature, /?', and the 
'true' one, (3. 

For the record, we give the summed forms of the scalar and spinor free energies 
on the two- sphere, 

F> (R\ = -— V cosh ("W 2 ) (oa) 
fscVV 2(3 ^ ims inh\m(3/2) 1 J 

F *p(P) = 'a E • \ 27 a /oV ( 4 °) 
P msmh (rap/2) 

We must now turn to the n-summation in (25). The calculation is much the 
same as for the spinor field. 

The prolate form obtained by the Watson-Sommerfeld method is, 



m = ZA1 C (3) - 2n r ! A ~ B 
sc " 8tt 2 ^ ] J exp(27rt) + 1 exp(27ryt)-l 



This is handy if we want the extreme prolate (large Z 3 ) limit. From the first term, 
the leading behaviour of the entire scalar effective action is 

^sc-^ic 2) ~-0C fl (3), / 3 -oo (42) 

showing a sort of duality. 

The oblate expression is found to be, 

w (2) _ f°° ^t 2 dt \ [\.„. n ^ ur _.^ [ C ^x smh(2n Bt) -Bsm(2nxt) 



exp(2-7rt) + 1 



' dy y coth(7ryt) — / dx- 
Jo 



cosh(2irBt) — cos(27rxt) 

(43) 

and the total scalar free energy, from (25), is 

W sc = C -^ + W^+W^. (44) 
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For prolate deformations, Wac is given by (41) and for oblate by (43). There 
is no Wsc^ for the prolate case and it equals (29) in the oblate. 

The corresponding 'free energy' in both cases is defined to be $ sc (/9) = —W sc /f3 
and from (44) we see that $ sc (/3) ~ 3> S c\P) so that from (38) this quantity has the 
same high temperature behaviour as on 5" 2 . This statement is expanded in section 
9. 

Again we note that, in the A — and C = limits, both (41) and (43) reduce 
to just the first term of (43). 

5. Spinors revisited. 

We now return to (13) for the spin-half (Pauli) free energy on the squashed 
three-sphere and again decompose the sum into even and odd pieces, 

9 oo 1 oo 9 oo (_-\\m 00 

d>(3) (/3 ) = _4y _L y ne -™ np'/4_^y <__LL y ne-mnP/4 (45) 

m=l n=2,4,... m=l n=l,3,... 

which presents a sort of twisted situation with thermal fermions being associated 
with scalar modes on S 2 and vice versa. 
Let us rewrite this expression as 

to be compared with (34). 

To obtain the leading behaviour we can make use of the fact that the functional 
relation (37) takes the twisted form, 

QM = Q' SP W) - 2Qs P (2/3). (47) 
The coefficients of the C 's in (35) and (36) will then be switched i.e. 

0?(3) r sp_ 4C«(3) 
2tt/3 3 " f3 3 



and 



then from (46) 



z?' f/^ 3C fl (3) 3^(3) . . 

F - (,3) ~ - &^r c » = -^a^ (49) 



Again, this follows also from the summed form (12). 

Once more we see that the leading behaviour on the extreme oblate three- 
sphere corresponds exactly to the high temperature form on the two-sphere, but 
not so the sub-leading terms. 
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6. Poles, residues and coefficients. 



In [9] we determined the analytic structure of the spinor zeta function, and 
thence the coefficients in the heat-kernel expansion. The same procedure can be 
pursued of course in the scalar case and just needs the complete pole structure of 
C sc ( s )- The intermediate Plana form (25) shows a pole at s = 3/2 coming from the 
Riemann ^-function. This is the 'volume' pole. It has a residue 4Z 3 which agrees 
with the general value, 

Cm 



(4n) d / 2 r(d/2-m) 

for the residue at s = d/2—m (m = 0, 1, . . . , ^ d/2) in a (/-dimensional manifold, Ai, 
if we recall the Weyl result, C = \M\. For our squashed three-sphere, \S 3 \ = 167r 2 / 3 . 

To obtain all the poles, we can proceed as described in section 2 the only 
differences being the change in sign and the bosonic factor. In the scalar case ( sc (s) 
has poles at s = 3/2 — m, m = 1, 2, . . . with residues 

R m = (2/ 3 ) 3 - 2m (2 1 - 2 - - l)r m 

in terms of the spinor residues (14). Hence the heat-kernel coefficients are 

C m = \S 3 \(2 1 - 2m -l)(ll-l) m ^. (51) 

ml 

Removing the volume factor, and taking the Z 3 — > limit yields Mullholland's 
results for the unit two-sphere scalar coefficients, 

lim = (-l) m (2 1 " 2m - 1) 
i 3 ^o 1 53 1 m \ 

The spin-half case is similarly treated. 

Shtykov and Vassilevich [14] have computed the coefficients on the deformed 
unit three-sphere for the minimal operator, —A, To get these coefficients one can 
simply multiply the above expansion by the expansion of exp(£/4) in the usual way. 
For example 

4 — Zi ,~o,R 



/^minimal I r<3 I 3 I r<3 

Cl - |6| ^2~- |6| 6 



as required. 
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8. An alternative summation and a puzzling coincidence. 



It is always a good idea to pursue different paths to any required quantity. 
Apart from inspiring confidence, or otherwise, in the final answer, it often reveals 
unexpected subleties. 

Instead of applying the Watson- Sommerfeld method to the n-summations in 
(6), (8) and (25), we can proceed as Okada does, [7], and employ the method first 
used by Minakshisundaram when discussing the ^-functions on spheres and later, 
more extensively, by Candelas and Weinberg [15] for the same purpose. The idea is 
to rewrite the summand using the Laplace transform 

in order to effect the n-summation. 

2 2 

For oblate l 3 < 1, A becomes imaginary (A = —C ) and the Bessel function 
becomes a modified one, I s -\/2- However there is an obstruction to the immediate 
application of a formula such as (52) in that, for t > 1/C, C will be larger than at 
least one n, violating the conditions of the identity. Our attitude in this section is to 
take the prolate expression as far as possible, and then to continue in The form 
of the eigenvalues shows that the singularities of the (^-function, as a function of Z|, 
all lie on the negative real axis and therefore should not prevent the continuation 
through Zjj = 1. 

Inserting (52) into the middle term of (6) allows the n-sum to be done and we 
find _ 

T(s) J sinh 2 (^/2)^ {AZ) J s-^ AZ > ^53) 
omitting the overall t-integral for the time being. At this point, following Candelas 
and Weinberg, the z integral is replaced by a contour one in the cut, complex z-plane 
(this could have been done earlier) to get for (53) 

Vk2 1/2 ~ s f , zsmBz 2s _ 2 . . u/2 _, T , , s lrtS 

— ? N x , : / dz = z 2s 2 (Az) 1/2 s J q _ 1/2 (Az) (54) 

(l + exp(2z7rs))l» Jc sinh 2 (^/2) V } sl/2V ' V ' 

where C can be taken as, say, z = x + iY (— oo < x < oo). The choice of the 
constant Y depends on any singularities possessed by the integrand. In this case, 
there are poles at z = 2pin (p = ±1, ±2, . . .) and so Y should be less than 2n. 

At this point we can again check the form of the residues (14) at the poles of 
f{s), which, in the present representation, occur at the zeros of 1 + exp(2z7rs), i.e. 
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at half odd integer s. For s > 1/2, the original form, (53), is convergent so the only 
possible poles are at s = 3/2 — m, m = 1,2, — For these values of s there is no cut 
in the z-plane and, also, the integrand is an odd function of z. Hence, we can add 
the expression for the reversed contour, —C (z = x — iY, x running from +oo to 
— oo) and divide by two. These two contours combine to give a clockwise contour 
around the pole at the origin and the integral is simply evaluated by residues. 
For (54) we obtain the value, putting back the t-integration, 

Jo 



r(3/2 - m)(m - 1)! J exp(2nt) - 1 

and we regain (14) which expression can be continued without ambiguity into the 
oblate region. Of course, being geometrical, the final result must be valid in both 
the prolate and oblate cases. 

Our main concern is with the derivative at zero, f'(0). In (54) one must 
differentiate the l/r(s) factor and set s — in the rest to obtain a non-zero answer. 
We find for the corresponding contribution to /'(0), 

f°° dt f dz — — 

/o(0) = - / 7^—, - / o sin Btz cos Atz (55) 

J2K 1 J exp(27rf)-l J c zsmh 2 (z/2) V ; 

and now the t-integral can be done (interchanging the limiting processes) using the 
standard formula 



I 



'o 

to produce 



" , SmCit 1 /«N 1 IT ^ 

dt : r r = - coth a/2 , Im a < 2%, 

exp(27rt) - 1) 4 v / ' 2a 1 1 



f(o) 1 / dz \ V' ( - 1 

2{ ] ~ AJc zsinh 2 (^/2) [ V ((B + A)z - 2ivip) + 



■p= — oo 



((B + A)z - 2nip) ((B - A)z - 2nip) 



If dz ( z sinh Bz 



4 Jc z 2 sinh (z/2) \ coshS^ — cosh Az , 

1 c (<\\ - 1 / dz ( sinhS^; \ 
n 2 flU ^ J c zsiYih 2 (z/2) VcoshS^ - coshA^/ 



27r 2 n 4 Jc z sinh 2 (z/2) \cosh Bz - cosh 

(56) 

The integrand is even in z ensuring that the integral is real. Convergence at large \z\ 
is also assured. We also note that the Cr(3) term cancels against the contribution 
to /'(0) from the first term in (22). 
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The price to be paid for performing the t-integration is a further constriction 
of the contour C, for which, in (54), Y lies between and 2rc. The condition on 
the validity of the t-integration means that now < Y < 2n/\B ± A\, or, if we set 
Z 3 = cosh 7, < Y < 7r(coth7 ± 1). 

The term in brackets in (56) presents poles at 

4±) = ^==ipn(l±^^=)=tpn(l±coth 1 ), p = ±l,±2,.... (57) 
B ± A a// 3 — 1 

Starting from the prolate side, all the plus sign poles, z p + \ p > 0, in (57) 
lie above the lowest existing pole at 2-rri and cause no problems. However the 
minus sign, p = — 1 pole, z_ x , falls below this pole when l 3 > 3/(2^/2) pa 1.06066 
and approaches the origin as / 3 increases, as do all the other minus sign poles, 
symmetrically with the sign of p. The contour C has to be adjusted to lie below 
this pole. This is what the above condition means. 

Just as a check, we have numerically evaluated (55) and (56). Typically, for 
7 = 0.5, (Z 3 = 1.1276) we find that (55) gives /£(0) « -0.123, for values, Y = 3 
and Y = 4, which straddle the pole at z = z7r(coth7 — 1) rs z'3.6565, the order of 
integration making no difference. In contrast, (56) gives —0.123 if Y is below this 
pole, but 0.143 if above. The difference is accounted for by the pole residue, as has 
been numerically checked. 

As the oblate case is approached, l 3 j 1, the plus sign poles, z p + \ with positive 
p run away up the imaginary axis to zoo, as do the minus sign poles, z p , with 
negative p. Conversely, the plus sign poles with negative p, together with the minus 
sign poles with positive p, run off to —zoo. 

As Z 3 passes into the oblate regime, the poles reappear at the complex positions 



(±) =pn(i±tan6) (58) 



z p 

We should follow them round carefully. It is apparent that some prolate poles with 
positive (negative) p, disappear with negative (positive) imaginary parts, only to 
reappear as oblate poles with imaginary parts of opposite sign. They have thus 
crossed the contour C (at complex infinity) and, if we want to leave C unchanged, 
it is necessary to include the corresponding residue contributions. These yield, 
choosing — 1 = i^l — /|, 



1 oo 

IE 



2 =1 psinh 2 z p 
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Symmetry under A — > —A (or C — > — C) is maintained automatically because 
sinh 2 Zp - " 1 = sinh 2 Zp + \ in tune with the fact that the branch point in (57) at / 2 = 1 
is purely artificial. 

The complete oblate contribution /^(O) is then 



m = ^c«(3) - \ j 



dz ( sinh Bz 

-=s=- ) + 



- 

oo 



c z sinh (2/2) V cosh Bz — cos 



1 



\ 2m sinh 2 (mar tan #) ^ (2m + 1) cosh 2 ((2m + l)7rtan(0)/2) 

(59) 

Other values of interest are /(n), with n a positive integer. These are needed 
in the expansion (5). Looking at (8) the first two terms are easily evaluated and 
the last one vanishes, as shown earlier. The remaining one, f2{s), is under present 
investigation and we find from (54), 

/a(n+ i) = r ' j\ r dz ^^ z ^(j tz) -n-^ j ^ 

n\ J exp(27rt) - 1 J smh (z/2) 1 

(60) 

This time the z integration can be done. As an example consider / 2 (1) 

= ^ .oc ^sin^)sin(^) 

2V ^ 7o exp(27rt) - 1 J Q Asinh 2 z/2 

and use the standard integral 

[°° , sin ax sin 6x „ 6sinh(2a7r) — asinh(267r) . . 

/ dx 5 = 2tt -) { /, / 61 

Jo sinli x/2 cosh(2ayr) - cosh(26yr) 



to leave a single numerical quadrature. In the oblate case (A — > iC), the integral 
in (61) always converges. 

The other values can be reduced similarly to involve integrals obtained by 
repeated differentiation of (61). 

For these values of s > 3/2, one can also calculate f(s) by direct summation 
which, purely numerically, is probably more convenient. 

For scalar fields, the evaluation of the n-summation term in (25) is precisely 
the same as decribed above for the spin-half case. The only changes are the overall 
sign of the term and the 'bosonic' factor l/(exp(27r£) + 1), which results in some 
algebraic differences. In place of (55) we have (remember W = — £'(0)/2 here) 

Wj® = ~[ / % sin Btz cos Atz (62) 

2 J Q exp(27rf) + 1 J c z sinh 2 (z/2) 1 ' 
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and this time we use the integral 



v , sinat 11 , . 
at — — r —r = — -cosech (a/2), |Ima|<27r, 



, exp(27rt) + 1) 2a 4 
to give in the prolate case, 



v y sc 



1 r dz sr^ 

~SJc ^sinh 2 (^/2) ±^ 



OO 

dz / (-l)P (-l) p 



c z sinh 2 (z/2) v ^ \ ((B + A)z - 2mp) ((B - A)z - 2nip) , _ , 

1 1 r dz sinh(S^/2) cosh(A^/2) 

4yr 2 flU 4 7 C 2sinh 2 (^/2) coshS^-coshA^ 

with the same restriction on the contour, C. 
The oblate expression is found to be 

dz sinh ( Bz/ 2 ) cos (Cz/ 2 ) 



c ^ sinh (z/2) cosh B z — cos C z 

y — ^ +y 



4 V ^ 2msinh 2 (m7rtan^) ^ (2m + 1) cosh 2 ((2m + l)7rtan(0)/2) 

(64) 

When this is substituted into the total expression (44) we again note the cancellation 
of the Cr(3) terms. 

There appears to be no ambiguity in the prolate calculation, and we have 
precise numerical agreement between (56) and (20) for spinors, and between (63) 
and (41) for scalars. 

However, extra pole terms seem to arise when the Bessel technique is 'contin- 
ued' from the prolate to the oblate regime, as exhibited in (59) and (64). These 
terms do not affect the perturbation expansion and indeed the first lines of (64) and 
(59) equal the whole of (43) and (19) respectively. Since they involve just a single 
quadrature, they at least provide better numerical alternatives. 

The puzzle is, partly, that these extra terms precisely cancel the oblate contri- 
bution, W S c (29) , in the scalar case, and double this up in the spinor case. This is 
a curious coincidence that we cannot explain. 

If one did believe the forms derived in this section, then there is no divergence 
in the total effective action as l s — > 0, for scalar fields. This would seem to run 
counter to a general expectation, derived from experience with the Selberg-Chowla 
formula, and also shows an unbelievable difference between scalars and spinors. 

It is clear, however, that the results on the C = case based on a perturbation 
expansion of the original sum form about C = (or about A = 0), discussed in the 
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earlier sections, show that there is no divergence in as l 3 — > and hence that 
the results presented in the previous sections of this paper are the proper ones. We 
therefore discard the second lines of (64) and (59). Doing this allows one to find 
the L — > limits of in closed form, 



and 



rr sp 



dz 1. „ 

3 = _ln2 

c z sinh z cosh z 4 



8tt 2 



C«(3) (65) 



Itc 8 J c z si 



cosh 2z dz 



1 n 1 
- In 2 H pt 

2 4tt 2 



Ca(3) (66) 



sinh 3 2; cosh z 

where the contour has been pushed upwards through the poles of the integrand and 
the residues summed. 



9. Preliminary comparison with the bulk results. 

We are now in a position to gather our results together and give the total high 
temperature forms of the free energy, defined by $ = — W//3. 
For the scalar field, from (44), (65) and (38), 

~ + ^ + \ (j^CP) - iln2 + C'(-l)) , (67) 

while, for the spinor field from (22), (66) and (50), 



*.p(« ~ - ± W - \ (^0,(8) + \ ln2 + 2C'(-1)) • 



(68) 



According to [1] and [2], the Taub-Bolt-AdS bulk, four-dimensional entropy 
goes like 1//3 2 at high temperatures. This behaviour agrees, of course, with our 
expressions for the free energy, or action, with which we prefer to work. 

Extracting the leading terms from [1] 

' \j3 2 8jt 2 1024tt 4 ' ) [ ' 

where a is a constant which does not concern us here. The identification parameter 
k has been included and should be set equal to unity to conform to the assumptions 
of the present paper. In this case, (69) agrees with the calculations in [2]. 

We notice immediately the absence of any logarithmic or transcendental terms. 
These can be eliminated from our expressions by choosing a combination of two 
scalars plus one spinor, yielding the behaviour 



X * (3) (^-l08V 



two scalars + one spinor 

which is the best that can be done with the fields available 



(70) 
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10. Conclusion. 

We have presented formulae for the determinants of spinor and scalar fields 
on the squashed three-sphere and have determined explicitly their leading 1/ f3 2 be- 
haviours which exactly correspond to those on S 2 x S 1 . The non- leading behaviours 
have also been found where the twisted nature of S 3 shows up. The scalar expres- 
sion is a mixture of thermal scalar and spinor on S 2 while the spinor form combines 
twisted scalar and spinor, i.e. fields with the 'wrong' thermal periodicity. 

The next step in the calculation, which will be exposed at a later date, is 
to place k identifications on the ^-circle. Since the resulting manifold is more or 
less locally unchanged, one would expect the highest terms to remain the same, 
depending, as they do, on the S 2 geometry. This is born out in (69) and (70), and 
is a comment made also in [2]. 

The extreme prolate limit presents certain technical difficulties which have still 
to be addressed. The action for the branch of the bulk action that leads to (69) 
tends to a constant as Z 3 — > oo while the other branch diverges like Zf, in agreement 
with (42) for the scalar field. For the spinor field, one needs to go beyond the 
expansion in Z 3 and to allow for sign changes in the spectrum. For the scalar field 
we should also compute for the minimal Laplacian. 

It should be noted that our W corresponds to —I of [1,2]. I is the difference 
between the actions of Taub-Bolt-AdS and Taub-Nut-AdS and is only defined for 
certain regions of Z 3 . There appears to be no corresponding limitation in the present 
calculation since the round case, l 3 = 1, is always accessible, whatever the operator. 

Finally we remark that it is necessary to consider vector fields, for which the 
eigenvalues are also available, [9,11]. 
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